Interplay between strain and oxygen vacancies in lanthanum aluminate 
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We evaluate the interplay between epitaxial strain and oxygen vacancy formation in the 
perovskite-structure oxide lanthanum aluminate, LaAlOa. Using density functional theory within 
the GGA+C7 approximation we calculate the dependence of the oxygen vacancy formation energy on 
the biaxial strain conditions. We find that the change in formation energy with strain is negligible 
over the range of strain values usually accessible through coherent epitaxial growth. Our findings 
suggest that experimental reports of different oxygen vacancy concentrations in strained films result 
from extrinsic factors, or from other impurities such as defect complexes. 



I. INTRODUCTION 

Perovskite-structure oxides, ABO3, are of fundamen- 
tal and potentially technological interest because of their 
numerous and often coupled functionalities, such as fer- 
roelectricity, magnetism, multiferroicity, metal-insulator 
transitions, etc. In recent years, it has been shown that 
functional behavior in such oxides can be enabled or en- 
hanced when they are prepared in thin-film form. For 
example, ferroelectricity has been reported in thin films 
of SrTiC>3, which is a quantum paraelectric in its bulk 
form, and enhanced polarization has been observed in 
thin films of ferroelectric BaTiC^^ In both cases, the 
improved behavior is attributed to the bi-axial epitaxial 
strain introduced when the film is grown coherently on a 
substrate with a different lattice constant. 

It is usually assumed that the change in in-plane lat- 
tice constant associated with epitaxial strain is accom- 
modated through structural distortions such as changing 
the internal bond lengths, or changing the angle or pat- 
tern of rotations and tilts of the oxygen octahedral— 
An additional possibility is that the film changes its de- 
fect profile, for example through changing the concen- 
tration of oxygen vacancies or the degree of cation non- 
stoichiometry. In particular the presence of oxygen va- 
cancies in perovskite oxides is believed to result in an ex- 
pansion of the crystal lattice due to underbonding caused 
by the two missing oxygen electrons . 27 ' 28 This suggests 
that the increase in volume associated with biaxial ten- 
sile strain might also increase the concentration of oxygen 
vacancies. Such changes are of profound importance as 
the presence of anti-site defects, interstitials, vacancies 
and/or defect complexes can have striking effects on the 
properties of perovskite oxides. 

Any correlation between epitaxial strain and intrin- 
sic defect profile is difficult to quantify experimentally. 
Strain is introduced by changing the substrate, which 
necessarily leads to other changes in growth conditions; 
the different growth conditions can then cause changes 
in the concentration and type of defects that might be 
unrelated to strain. In this work we use first-principles 



density functional theory to investigate the relationship 
between strain and the formation energy of oxygen vacan- 
cies in the model perovskite-structure oxide, lanthanum 
aluminate, LaAlC>3 (LAO). Bulk, strain- free LAO has 
Ric symmetry with a slight distortion of the rhombohc- 
dral angle away from the ideal 60°. The oxygen octa- 
hedra rotate by ~6.3° in alternating directions around 
the pseudocubic [111] direction, in the Glazer 
notation tilt pattern»i2r— It has a large band gap of 5.6 
eVfi^ and so is robustly insulating even within the local 
density or generalized gradient approximations to den- 
sity functional theory, which are known to systematically 
underestimate electronic band gaps. 8 This is particularly 
important for first-principles studies of defective systems, 
so that defect-related transition levels which are experi- 
mentally found to be in the band gap do not erroneously 
fall outside of the band gap, changing the predicted elec- 
trical characteristics. In addition, LAO is structurally 
robust with no tendency to Jahn- Teller or ferroelectric 
distortion. 

First-principles calculations were used recently to de- 
termine the intrinsic structural distortions induced by 
changes in the in-plane lattice constants of LaA103. 14 
It was found that the [lll]-oriented rotation pattern of 
bulk LAO is stable only for a small in-plane biaxial strain 
region (-0.2 to 0.1 %) and larger strains result in differ- 
ent patterns. Biaxial compression changes the rotation 
scheme to a°a°c~, in which the octahedra rotate alter- 
nately along the [001] direction, and tilts around the in- 
plane axes are suppressed. For biaxial tension, the oc- 
tahedral rotation scheme changes to a _ a _ c°, with alter- 
nating rotations around the [110] direction. 

First-principles calculations have also been used to ex- 
plore the structure and formation energies of various 
defects in the high temperature, high symmetry cubic 
Praim structure of LAOi 15 ' 16 Luo et al. calculated the 
atomic structures, energy levels, and formation energies 
of charged interstitials, vacancies, Frenkel pairs, antisite 
defects, and Schottky defects. As expected, they found 
that the formation energies of O vacancies are lowest in 
reducing conditions, whereas the La vacancy formation 
energy decreases with increasing O chemical potential. 



They also found that considerable structural distortion 
occurred around many of the defects, with substantial 
changes in the inter-ionic distances. It is unclear from 
Ref. [l6| whether the lattice parameters were fixed to the 
ideal values during the calculations; changes in cell pa- 
rameters associated with the different defect types were 
not reported. Likewise the focus of Xiong et al. in Ref.[l5l 
was on the electronic energy levels generated by the for- 
mation of defects, rather than on any associated changes 
in lattice parameter. 

Here we present results of our first-principles study 
combining both epitaxial strain and defects in our calcu- 
lations. Our goal is to determine whether a critical value 
of strain exists at which it is intrinsically more energet- 
ically favorable to incorporate defects into the sample, 
instead of further deforming the bond lengths and angles 
of the ideal stoichiometric structure. 



II. COMPUTATIONAL DETAILS 
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FIG. 1. Calculated oxygen octahedral rotation angles within 
the LDA+U and GGA+U methods as a function of the Hub- 
bard U parameter. The experimental value is shown by the 
horizontal dashed line. 



Our calculations were performed using density func- 
tional theory as implemented in the Vienna ab initio 
simulation package (VASP).— The projector augmented 
wave (PAW) method was used for core-valence separa- 
tion with the supplied VASP potentials for the La, Al 
and O atoms' Test calculations for bulk LaA103 were 
performed using a 10 atom rhombohcdral cell (pseudocu- 
bic angle fixed at 90 degrees) with a 4x4x4 Monkhorst- 
Pack fc-point mesh and a plane-wave basis energy cutoff 
of 600 eV. We defined biaxial strain as 



where ag is the calculated equilibrium lattice constant. 
Once strained, the out-of plane (c) lattice parameter was 
relaxed and the ionic positions optimized until the forces 
were below lmeV/A. 

A. Choice of exchange-correlation functional 

We first tested the behavior of ideal bulk LaA103 to se- 
lect an appropriate exchange-correlation functional that 
adequately describes both the structural and electronic 
properties. While the local density approximation (LDA) 
has been previously shown to give a good description of 
the structural properties of LAO^ it incorrectly places 
the highly localized La / states at the bottom of the con- 
duction band^ where they will likely spuriously inter- 
act with oxygen vacancy states in our defect calculations. 
One apparent solution is to apply a Hubbard U correc- 
tion within the LDA+J7 formalism; a U value of 11 eV 
was used previously to correct the La / state position in 
studies for which detailed structural distortions were not 
important^ We found, however, that, while the LDA ob- 
tains almost perfect agreement with experiment for the 
magnitude of the octahedral rotations, the addition of 
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FIG. 2. Calculated GGA+U density of states for LaA10 3 for 
different U. In all cases, the top of the valence band is set to 
OeV. As U is increased, the La / states increase in energy and 
shift away from the conduction band edge (red arrow). The 
valence band density of states is effectively independent of U. 



a U to the La / states strongly suppresses them giving 
unacceptable structural properties. In contrast, the gen- 
eralized gradient approximation (GGA), over-estimates 
the octahedral rotations and addition of a U to the La 
/ states reduces them towards the experimental values, 
while also correcting the energy of the La / states (Fig. 
[21) ■ Our tests of electronic properties and structural prop- 
erties as a function of U led us to a choice of U = 10.32 eV 
on the La / orbitals as the optimal choice, within PBE 
GGA 23 and the Dudarev method. 24 Double-counting cor- 
rections were made in the fully localized limit. 
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FIG. 3. The two symmetry-inequivalent positions for an oxy- 
gen vacancy in biaxially strained LaA103 within our supercell. 
In this study we use position 1, in which the oxygen vacancy 
lies in the normal axis to the strain rather than in the basal 
plane. 



B. Oxygen vacancies 

Oxygen vacancy concentrations of ^4% (one out of 
24 atoms) were studied using supercells containing 40 
atoms (a 2x2x2 supercell of the 10-atom rhombohedral 
unit cell with the pseudocubic angle fixed to 90 degrees); 
this allowed us to reduce our Monkhorst-Pack mesh to 
lxlxl. When biaxial strain is included, two distinct 
placements of the oxygen vacancy - within the biaxially 
constrained direction or in the unconstrained relaxed c 
direction. (Fig. [3]) - are possible within our choice of su- 
percell; here we studied the perpendicular configuration. 

The oxygen formation energy was also calculated for 
a larger 4x4x4 supercell to assess proximity to the dilute 
limit. The change in formation energy between the 2x2x2 
cell and the 4x4x4 cell (corresponding to a change in 
oxygen vacancy concentration from ^4% to ^0.5) was 
approximately lOmeV. 

III. RESULTS 

A. Calculated structure of unstrained LaAlOa with 
oxygen vacancies 

First we introduce an oxygen vacancy into the 40 atom 
supercell and fully relax the atomic positions with no 
strain applied, constraining the lattice parameters to 
their bulk values. We find that the octahedral rota- 
tion pattern changes from the a~a~a~ pattern of bulk 
LaAlC>3 to the a~a~c°, in which the octahedral rota- 
tions around the c axis are lost. Interestingly this is the 
same pattern as found in defect-free LAO under biaxial 
tension. The deviation from the bulk rotation pattern 
with the introduction of the oxygen vacancy is likely re- 
lated to the long range coherence of octahedral rotational 
instabilities, and their interruption by the presence of 
the vacancy. This finding that an oxygen vacancy has 
long range structural implications opens the possibility 
for strain to change the formation energy for oxygen va- 
cancies. 
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FIG. 4. Calculated densities of states for unstrained LAO 
with (black line) and without (red line) oxygen vacancies. 
The densities of states are aligned so that the tops of the bulk 
valence bands lie at eV. The occupied states are shown with 
shading; in the vacancy structure, the impurity states are fully 
occupied and lie within the band gap, 2.6 eV from the valence 
band edge. 

B. Electronic properties with oxygen vacancies 

Our calculated electronic densities of states for un- 
strained LAO in both ideal and defective structures are 
shown in Figure |4j The oxygen vacancy state is fully 
occupied and falls within the band gap of the ideal bulk 
material. Apart from this isolated impurity band, the 
ideal and defective densities of states are very similar. 
The different k-point sampling between the two calcula- 
tions causes the slight change in the smoothness of the 
plots. 

C. Combined effect of strain and oxygen vacancies 
on structure 

We find that the introduction of an oxygen vacancy 
into the strained structures does not qualitatively change 
the octahedral tilt patterns from those of the corrsepond- 
ing strained structures without defects. Under biaxial 
compression, the defect structure adopts the a°a°c~ tilt 
pattern found in the defect-free structure under biax- 
ial compression; under biaxial tension the a~a~c° struc- 
ture of the unstrained oxygen vacancy structure, and the 
defect-free tensile-strained structure is retained. 



D. Strain dependence of oxygen vacancy formation 
energy 

Next we calculate the formation energy of oxygen va- 
cancies as a function of biaxial strain. The formation 
energy, fi, of an oxygen vacancy can be calculated from 

n = E tot (VS)-E tot (bulk)+Nno (2) 
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FIG. 5. Calculated oxygen vacancy formation energies, Q, 
as a function of oxygen chemical potential, /io for a range 
of strain states. Lower chemical potentials correspond to 
oxygen-poor conditions and vice versa. Results for four com- 
pressive strain and four tensile strain values (± 0.5%, 1%, 
1.5%, and 2%) are shown. The bulk (unstrained value) falls 
lies with the tensile strain values. 



where E t ot{Vo) is the energy of the supercell with an 
oxygen vacancy, E to t(bulk) is the energy of the defect-free 
supercell, N is the number of oxygen vacancies (one in 
our case) , and is the chemical potential of the oxygen 

• 2Q 

reservoir— . 

The oxygen chemical potential (po) is defined as 

Mo = ^OMuik - ^E(0 2 ) (3) 

where [io,buik is the bulk chemical potential and E(Oi) 
is the formation energy of an oxygen molecule. It is clear 
that the formation energy depends linearly on the oxy- 
gen chemical potential; at low oxygen over-pressure, the 
formation energy is lower yielding a higher concentration 
of oxygen vacancies and vice versa. 

In Fig. [5] we show our calculated formation energies at 
different biaxial strain values (from -2% to +2%) as a 
function of the oxygen chemical potential. Note that the 
formation energies are large - > 7eV - consistent with 
previous results for related oxides . 15 ' 16 As expected from 
our earlier simple volume arguments, the formation ener- 
gies at a given chemical potential are higher under biaxial 
compression than under tension. The strain dependence 
of the formation energies is small in comparison to the 
total formation energies, with a 2% increase in compres- 
sive strain increasing the formation energies by only ^20 
meV. 



We find that the strain dependence of the formation 
energy at a specific oxygen chemical potential is much 
stronger for compressive than for tensile strains. This 
is likely a consequence of our choice of location of the 
oxygen vacancy within the unit cell: Under compression 
rotations around the [001] direction dominate the struc- 
ture, and these are more strongly affected by the in-plane 
oxygen vacancy. 

IV. DISCUSSION 

In summary, we have shown that density functional 
theory within the GGA+U approximation gives a good 
description of both the structural and electronic proper- 
ties of LaA103, when a U of 10.32 eV is applied to the La 
4/ states. Perhaps surprisingly, introduction of oxygen 
vacancies at 4% concentration in axial positions induces 
a change in the tilt pattern from the bulk a~a~a _ to the 
a~a~c° pattern; this is the pattern that was previously 
identified for defect-free LaA103 under tensile strainJ^ 

We found that the formation energy for an oxygen va- 
cancy in LaAlC>3 is between 3-10 eV depending on the 
oxygen chemical potential. The change in this value with 
experimentally accessible strain values are orders of mag- 
nitude smaller, with a ^75 meV increase for 2% compres- 
sive strain and 5 meV decrease for 2% tensile strain. In- 
terestingly, the change in vacancy formation energy with 
strain is comparable to the change in intrinsic total en- 
ergy (Fig. 2 of Ref. Il4l ). suggesting that under suitable 
conditions strain could be partially accommodated by 
change in vacancy concentration. It is also likely, how- 
ever, that strain will modify barriers for defect migration, 
making a direct experimental comparison of equilibrium 
defect concentrations at different strain values fundamen- 
tally prohibitive. 

Finally we mention that in perovskite-structure (and 
other) oxides in which the B site is a transition metal, 
the introduction of oxygen vacancies causes a change in 
oxidation state of the B site. Such changes in oxidation 
states are accompanied by large changes in coordination 
volume. As a result, stronger strain-dependence of oxy- 
gen vacancy formation energies should be expected. 
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